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ABSTRACT 


A three dimensional finite element code is written for 
truss analysis and design. Trusses may be designed for mini- 
mum weight subject to constraints on: member stresses, Euler 
buckling, joint displacements and system natural frequency. 
The optimum configuration may be found in addition to optimi- 
Zation with respect to member sizes. 

The finite element code may be used as a stand alone 
analysis tool or may be coupled to an optimizer of the user's 
choice. The finite element displacement method of analysis 
is used for static analysis and eigenvalues are calculated 
using the subspace iteration technique. 

A major portion of this work is the software (user guide) 
which 1s presented in detail with examples and results. Ex- 
planation of how this code is coupled to an optimizer is 
given. 

Design variables are member area sizes, joint coordinates, 


or both. Examples are presented to demonstrate the method. 
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Pa ODUCT TON 


Design optimization requires the minimization or maximi- 
zation of some parameter. Optimization of structures has had 
continuing changes since its development in the early 1960's 
with an active area of research being elastic truss structures. 
The main goal is to design structural systems that efficiently 
perform specified purposes. For any design to be acceptable, 
it must satisfy a variety of physical, aesthetic and economic 
constraints. 

Since most physical problems can be modeled by some mathe- 
matical formulation a computer program can be written to per- 
form the necessary calculations. 

The purpose of this research was to develop a FINITE 
ELEMENT code that could be easily coupled to an optimizer, 
thus (1) allowing testing and comparison of various optimizers 
and (2) provide a useful design program in its own right. 

The design problem considered in this study is the opti- 
mization of three-dimensional indeterminate trusses, for 
multiple static load conditions. The objective is to mini- 
mize the weight of the structure where the design variables 
are member sizes and joint coordinates. Constraints include 
stress, displacement, Euler buckling, and the natural fre- 


Guemey Of the system. 
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This document describes the use and capabilities of the 
finite element computer code to be coupled to an optimizer. 
The user's manual presented in Chapter V contains a simple 
design example in which the program is coupled to the CONMIN 
optimization code [Ref. 1]. 

Additionally, guidelines for coupling the code to an 
optimizer of the user's choice are presented. 

Several examples demonstrating the program under a variety 
of conditions are presented. Conclusions and recommendations 


for future work are given. 


i 





wie Oe IMT ZAt TON 


A. INTRODUCTION 

The main goal of structural engineering optimization is 
to design structural systems that efficiently perform speci- 
fied purposes. Selection of a specific algorithm must include 
the following considerations: 1) the structure should be 
analyzed as few times as possible, 2) the algorithm should 
minimize specific gradient information required, and 3) it 
Should provide reasonable assurance that an optimum design 
will be reached. 

The next few sections discuss the general formulation of 


the algorithm with respect to the above requirements. 


B. FORMULATION 
Minimum weight design of trusses is presented in the 
general form of a mathematical programming problem as follows: 
Masqamaze F CX) GaGa nk) 


Sulbspiect to: 


where F(X) is the objective function, in this cas@, weight of 
the structure be minimized. X is the vector of independent 
design variables and here contains the member cross-sectional 


areas as well as the coordinates of the joints. The inequality 
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constraints, Bal feo, mMvst be satisfied for the design 
to be accepted as feasible. These include limits on stress, 
Euler buckling, joint displacements and the first fundamental 
frequency of the structure. Side constraints yi and Xe are 
lower and upper bounds on the design variable. These may be 


treated as general inequality constraints. 


eee PARAMETERS 

The design variables, constraints, and objective function 
considered in the design process are discussed here. 

1. Design Variables 

The design variables the user has available are member 

cross-sectional area, reciprocal of member cross-sectional 
area, joint coordinates, or both member areas (or their reci- 
procals) and coordinates. The cross-sectional areas are AL 
k=1,NE, where NE is the number of elements in the truss 
structure. The joint coordinates are Set ieee ja Ne 
where 1 is the coordinate axis, j 1s the joint number and NJ 
is the total number of joints. Design variable linking is 
allowed in both member sizing and coordinate design variables. 

2. Objective Function 


The objective function considered here is weight. 


NE AL 
Weight (W) == p. "4 (Eq. 4) 


where 9p is the material density, A is the cross-sectional area 
and L; is the length of the member. The truss may be made up 
of members of differing materials. 


1a 





3. Constraints 





Th, Ses Chan. 5) 
Cae 
Stress: —+ -1< 0 Ta, NLC 
Co 
Sil we Ceamaro 
om 
nl Bie 
=i -1< 0 jesh NILE 
Oo 
where 0 1s the stress in member 1 under load condition j. NE 


is the number of elements and NLC is the number of loading 


Gemditions. 
¢ is the lower bound on the stress (maximum compressive 
oe ae 
stress) ando is the upper bound on the stress. The upper 


and lower bounds on stress may be different for each member, 
but are taken to be the same for every loading condition. 
Pe eee ets ha ie. 
The stress at which Euler buckling occurs is given by: 
a ee St The 1 NS (Eq. 7) 
7 2 


ley 
aL 


where the subscript 1 corresponds to the member number, Ee 
is Young's modulus and K. is a constant depending on the 
cross-sectional geometry of the member. 
5. Displacement 
Displacement limits are imposed at prescribed joints 


to create a constraint equation as follows: 
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mms 2 1 <0 Seg 8 
Mick 
ise 

Be apie 
—J>-1< 0 Gan?) 
“aig 


where Us sx and Us ax are lower and upper bounds on the displace- 
ment at joint 1 in the coordinate direction j under loading 
@emaition k. 
6. Frequency 
The first fundamental frequency of the structure is 


required to exceed the specified lower bound, so that 
imeem Xk <0 (Eq. 10) 


Poor oneeome ets ameexcelient source for the basic 


structural design formulation. 
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Peete ELEMENT METHOD 


A. INTRODUCTION 

Several features are desirable when the finite element 
methods of analysis are used in design optimization. First 
mremmumoer Of analyses for the structure should be kept toa 
minimum. Second, the amount of gradient information required 
during the design process should be reduced to shorten run 
times and computer storage requirements. Third, the user 
should be able to specify only that gradient information 


desired. 


Bee ANALYSIS 
Initial problem formulation includes member sizing, mate- 
rial properties (which may be different for each member), a 
set or sets of external loads, and specified support condition. 
The analysis for the stresses and deflections must satisfy 
the conditions of equilibrium of forces at the nodes and com- 
patibility of deformation. In this analysis the weight of the 
individual members are not part of the specified load conditions. 
Additionally for this analysis the following assumptions 
are made. Trusses will be treated as discrete elements, and 
each element will be treated as pin-connected with loads and 


reactions supported at the joints. 
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The general method of solution is as follows. The Dis- 
placement (Stiffness) method reference 2 considers the joint 
displacement components as the unknowns, and written in 


matrix notation for the general case is 


where K Peo eet sett iness Matrix, P is the vector opr 
@eerors Or applied ieads, and u is the vector of vectors of 
displacements. 

Once displacements at every node are known, the internal 
Bepees and stresses are calculated by applying the appropriate 
force-deflection relationships. 

When the system's natural frequency constraints are con- 
Sidered the design process requires the solution of an elgen- 
problem. This solution will determine the natural frequencies 
and normal modes of the structure. For linear elastic struc- 
tures, the finite element approach leads to the following 


equation of motion with free vibration conditions, 
Mor + Kaie= 0 Clan e2) 


where M is the global mass matrix, and u is the linear accel- 
eration vector. This leads to an eigenvalue problem of the form 
1K-\ 1d=0, and is solved here by the subspace iteration method 


reference 2. 


Ly 
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Seen pieNtTsS OF CONSTRAINTS 

Gradient computation of necessary functions in a design 
Optimization process arise from the need for derivative infor- 
mation for efficient mathematical programming. For structural 
optimization the gradient of displacement with respect to the 
design variables, du/dx, 1s needed from which VG, is calculated 
for stress, displacement, and buckling. 

Pemseder the finite element method Ku = P. Taking the 
derivative of both sides with respect to the design variable 


Xqg, we have: 


— [Ku] = — (P) Cedar 


Assuming the loads P are not a function of Xo 5 


MG) Ga) + kK (au) = 0 (Eq. 4) 
9X = = IXo 
; ) ~1 ) , 
Finally, we arrive at —~— (u) = -K [——- (K)] u (Eq. 5) 
IX, = 9Xo = . 
where Tee is the inverse of K. For efficiency an 1S NOs 


actually calculated. 
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TV. PROGRAM FEATURES 


A. INTRODUCTION 

Each computer code has its own special features with which 
the user should be familiar if the program is to be used 
effectively. The next few sections will discuss an overview 
of the code and a typical problem that might be solved using 
the program. This problem along with other numerical examples 
will be presented in detail with results in Chapter VI. 

The FINITE ELEMENT code was written to be used as a stand 
alone analysis program or an analysis code that could easily 
be coupled to an optimizer (of the user's choice) through 
Simple modifications to the main driver program. 

With user supplied area and coordinates coupled with input 
control parameters, the analysis mode will calculate the weight 
of the truss structure. Design variables may be chosen as 
member areaS, joint coordinates, or both. Additionally, 
gradient information will be calculated with respect to area 
variables, coordinate variables or both. Coupled to an opti- 
mizer the code will optimize the weight of the structure and 
print the final optimization information and gradient vectors 
of specified constraints. 

The following example of a 25-bar space tower presented 


in Tables (I-IV) shows some of the options the code contains. 
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fable III 


29-BAR TOWER SAMPLE OUTPUT 


AREA/COORDINATE VARIABLES 


THERE ARE 


per, at VA 
e20000 


AREA DESIGN VAR TABLES 


O-20000EtOL O-Z20000E+O01 O.20000E+OL O.20000E+01L 0 .20000E+01 


04 0-LODD0E-0% 0.10000E-04 0.-10000E-04 0.1L0000E-04% 0.LOD00E-04 


UPPER BOUNDS 
Qe50000E+02 O-S5ODOD0EtO2 O0-50000E+02 0.50000E+t02 0.50000E+02 0.50000E+02 
THERE ARE 


vai VALU 
0.37590E+ 


5 COORDINATE VARIABLES 
ES 
+02 0.37500E+02 0-1LO000E+03 0. LO000E+03 0.1L0000E +03 


LOWER BOUNDS 
O-LOOOCE+00 O-LOODDOE+00 O-LO000E+00 O.LO000E+00 0.10000E+00 
03 0.50000E+03 0.50000E+03 0.50000E+03 0.50000E+03 


JOINT OFSPLACEMENT CONSTRAINTS 
DIRECTION L=X_, 2=Y, 3=2, O=RESULTANT 


LOAD LOWER UPPER 
NODE OIR. COND. BOUND BOUND 
1 1 L -0.35)00E+00 Q0.350CE+00 
2 1 L ~0.3500E+00 0.350CE+00 
1 2 L -9.3500E+00 0.3500E+00 
2 2 1 -0.3500E+00 0.3500EF+00 
FREQUENCY CONSTRAINTS 
LOWER BOUNDS IN CPS 
FREQUENCY LOWER i 
NUMBER BOUND : 
L 0.1600E+02 
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POEs 


60 
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GRADIENT CF OBJECTIVE WAS CALCULATED 
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Ure UoeR GUIDE i 


A. INTRODUCTION 

In developing any computer code for engineering analysis, 
it is necessary to additionally develop concise, easily under- 
Stood software. This USER GUIDE is writen to be easily 
followed assuming minimal FORTRAN knowledge. The format 
moeerows that of the optimization code, COPES/CONMIN, 
reference l. 

This chapter is devoted to acquainting the user with the 
code and necessary input data. A simple 3-bar truss analysis 


1s used as the example. 


B. GENERAL FORMULATION 


Each data card is set up to show the input data necessary 


with the 3-bar truss example underneath. Formats are of "I5" 
mee FaO.0" type. "I" formats must be right justified, and 
"Pr" formats must have the decimal point. The number of cards 


read per data block is listed at the bottom of each block. 


©. CONSTRAINTS 

Constraints are calculated and stored in the G vector as 
listed in the following chart. The total number of constraints 
NCON = 2*NDSPLC + NFREQ + NE*NLC. When any of the constraints 


are missing from the G vector, all constraints are moved up. 
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For exemple, if there is no frequency constraint, then a dis- 
placement constraint would fill the first location of the G 


mMeccor. 


NFREQ 


Ga ROeURNCY “CONSTRAINTS ) 


2 NiepEe % NC 


(Dis PHAGEMENT CONSTRAINTS) 


G a ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee 
MEMBER #1 STRESS AND 
Rear. BUCKLING 
MEMBER #2 STRESS AND 
LC #1 BUCKLING 
MEMEER ANE Sa neos AND 
Ie #NLC BUCKLING 





zZ9 





Ds EXAMPLE 

Mjestitelalstavyout Of the 3-bar truss is shown in Figure 
6.1. Stress constraints were imposed as well as constraints 
on Euler buckling, displacement, and first fundamental 
frequency. 

With geometry specified as per the figure, 2 independent 
load conditions (P1,P2) were imposed with 3 member sizing 
variables (Al, A2, A3) linked so that Al=X(1), A2=X(2), 
A3=X(1). 


Pees sepemeses/ Conditions 
The material used had a density (p)=0.1 lb/in’. 


Young's Modulus was selected as Oe psi. 

LGADS : Pie 2 = 0.000 (ibs. 

Mita oLEMeTS9-15,000 < G- < 920,000 psi i=1,3 

Symmetry was maintained by linking of variables so 
Al=A3. Initial design began with Al=A2=A3= 1.0 in. One non- 
structural mass of 500 lbs was attached at joint 2. Displace- 
ment constraints were imposed between -.02 and .02 inches in 
the x-direction with the first natural frequency limited to 
a value if Zon HZ . 

ep ue control Parameters 

The following input control parameters are given for 
ease of the following example. NE=3, NJ=4, NCJ=4, NMT=1, NLC=2, 
NEIG=2, NEIGl=2, NFMASS=1, NEUBC=1, NDSPLC=2, NFREQ=1, LMASS=0, 
mp ebe-2 Table VY is a listing of commonly used nomenclature. 


é 
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TABLE V 


COMMON VARIABLE NOMENCLATURE 


A- MEMBER'S CROSS-SECTIONAL AREA 
BL- LOWER BOUND ON DISPLACEMENTS 
BU- UPPER BOUND ON DISPLACEMENTS 
DIR- DIRECTION 1=X%, 2=Y, 322 

O=Limit of RESULTANT DISPL. 
E- YOUNGS MODULUS 
ELNO- ELEMENT NUMBER 


EPSEIG- CONVERGENCE TOLERANCE OF EIGENVALUE 
OTIONS (DEFAULT=.000 ) 
FC1,FC2,-...FCN- LOWER BOUND ON FIRST,SECOND, ETC. 
SYSTEM NATURAL FREQUENCY 
FX- LOAD PORCES BEING APPLIED IN THE X DIRECTION 
FY- LOAD FORCES BEING APPLIED IN THE Y DIRECTION 
PZ- LOAD FORCES BEING APPLIED IN T DIRECTION 


HE 2 
GRAV- ACCELERATION DUE TO GRAVITY (DEFAULT= 386.4) 
ie CONSERALNT EDENTIFIER. IF NON-ZERO THE X-DOF 
IS CONSTRAINED 
fee weeNorRaLNT IDENTIFIER. IF NON-ZERO THE Y-DOF 
me CONSTRAINED 
ime GONSLTRAINT IDENTIFICR. IF NON-ZERO THE Z-DOF 
IS CONSTRAINED | 
JN- JOINT NUMBER (GLOBAL) 
REULER- EULER BUCKLING COEFFICIENT FOR BAR 
ELEMENTS 
Ee= LOAD CONDITION 
PRASS- LUMPED MASS OPTIONS. IF LMASS .NE. 0 LUMPED 
MASS MATRIX ED 
ims LWASS =O CONSISTENT MASS MATRIX USED 
NEIGI-NOMBER OF EIGENVALUE/SEIGENVECTOR TO BE 
EVALUATED DEFAULT eee aes NEIG+8) 
NEUBC- BUCKLING CONSTRAINT IDENTIFIERS . IF NEUBC. 
NE.Q0 EULER BUCKLING WILL BE IMPOSED IN 
BAR ELEMENTS 
NFMASS-NUMBER OF FIXED MASSES ATTACHED TO THE 
STRUCTURE 
NPREQ-~NUMBER OF FREQUENCY CONSTRAINTS 
NID- NUMBER OF INDEPENDENT DEGREES OF FREEDOM 


G 
” 


NLJ~ NUMBER OF LOADED JOINTS FOR THIS LOAD CONDITION 
NHT- NUMBER OF SEPARATE MATERIAL TY 

RHC- MATERIAL DENSITY 

SIGMIN~ MINIMUM ALLOWABLE ae 

SIGMAX- MAXIMUM ALLOWABLE ESS 

XA- INITIAL VALUE OF AREA DESTCN VARIABLE 

XAL- LOWER BOUND 

AU- UPPER BOUND 

XC- INITIAL VALUE OF COORDINATE DESIGN VARIABLE 
XCL- LOWER BOUNDS 

ACU- UPPER BOUNDS 
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The following user's manual is divided into blocks A 
through M. Appearing directly below each data field line 
are the parameters for the 3-bar truss example. It is impor- 
tant to note the user may choose any units; however, all units 


must remain consistent throughout the problem. 


rms) 














| 
| | 
| | 
| 

| | 
| | 
| 
| | 
| | 
| | 
| | 
pn nn nr rr | 
ir Ler | FORMAT | | 

| 20A4 | | 
Sn ee ee 

| -——-—_--_____— — = = | 
| 3-BAR TRUSS (EXAMPLE) | | 
j t— — ee | 
| 
| | 
| 
| | 
| | 
| - 
| | 
| 

| FIELD CONTENTS | 
| | 
| 

1 ANY 80 CHARACTER TITLE i 

MAY BE GIVEN ON THIS LINE | 

| | 
| | 
| | 
| | 
| | 
| | 
| | 
| | 
| | 
| | 
| | 
{ _ | 


2 


Th. 





mE ee oe [SS a See 


DESCRIPTION: Control Parameters 


mie ED EE a cee OS ee 


Format and Example 





“NEI NI} Neg] N ut] NL NLC FE zuBc] Nps?t Lc] Nar | [FORMAT - 
| { | Jers | 
| { { { | { 
u a =... 
|©6—Cl LL BRB 

[LL ee ae a a Sen = SS | 

{ NEIGINEIG1 ea NFREQ II vcLC] IRECP IRECP JForwat : 
mo | 2 1 | 1 2 0 | 515 

_ 

| | | | | | 
EE 


NOTE: DEFINITIONS OF INPUT CONTROL PARAMETERS FOR 
PROGRAM NEXT PAGE. 


: 

| 

| 

! 
ch 

iF 


| 
| 
| 
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NE-number of elements 
NJ-number of joints 
NCJ-number of constrained joints 
NMT-number of seperate material types 
NLC-number of load conditions 
NEUBC-buckling constraint identifier 
(If NEUBC.NE.0 -~EULER buckling con- 
straints will be imposed on bar elé.) 
NDSPLC-number of displacement constraints 
NER inewe orint control (if NPR1.ne.0 
input information will not be printed) 
NEIG-number of precise? eigenvalues to 


be evaluated 


NEIGI-number of eigenvalues to be evaluated 


DEPAULT=min. of (2"®NEIG , NNEIG#8) 
NFMASS~-number of fixed masses attached to 
structure 
NFREQ=-number of frequency constraints 
IDVCLC-design variable control parameter 
If (IDVCLC.EQ.1) NDV=NDVAR1 
If (I DVCLC.EQ.2) NDV=NDVAR1+ NDVAR2 
If (IDVCLC.EQ.3) NDV=NDVAR2 


ITRECP-recripocal variable identifier 


If (IRECP.GT.0) the X-vector contains 


the recripocal of area and VLB & VUB 
are changed accordingly 


SE se ee ee ee ee ee SS, ee OS GR: eM _ eg Ce Lee co es se CO: ce cs oe i ee ee 


| 
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ee a eee 


7 





the Lumped mass matrix is used. 


2 GRAV-accleration due to gravity 
(defaul+=386.4 inches/sec ) 


Laas | GRAV | EESEIG] | FORMAT 2 
| | 1) ae | 
| | | | 2F10.0 | 
es Sa 
0 | 386.4 | | | 
| : { | 
ee SS ae ae p | | 
| 
f 
i 
i 
| 
~ i 
| 
| 
| 
FIELD CONTENTS | 
1 LMASS-lumped mass option (if LMASS.NE.0) | 
‘ 
' 
{ 
( 


3 EPSEIG-convergence tolerance on eigenvalue 
solution. (default=.0001) 
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DESCRIPTION: Joint Coordinates §& 


Design Variable Linking 


Format and Example 


a nay 
|| San (oe oe il 0 | 0] 72010. 
haat 


JN- joint coordinate number 
X-x coordinate 

Y-y coordinate 

Z-z coordinate 

IX-design variable associated 
TY¥-design variable associated 
IZ-design variable associated 
PCX-participaticn coefficient 


Oo mon HD WN F&F WW NHN = 


PCY-partici pation coefficient 


2 
© 


PCZ-participation coefficient 


NOTE: Number of cards read=NJ 
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nl x Y | Z IX ty ]1z] PCX pce PCZ | Format 


|i 
| aisel 0 


em See 


. 
OTT ee ED a OL ee ETT OT A oi a AE so I ee i ee CO ES ne ee ee GEE ee eee a oo, eS coon, GE ete os ee eee ate ein: ae eee) eee SS ee ee a me ce ee eee a et 


CE 


Woche x eGOoOr dd’. 
wit hp y coomd: 
wit he ZoCcoomd:: 
Of  XK—-COorda. 
Of Ve'coord. 


of z-coord. 





SP AS ES ET A ED AALS ES APTS OE EES, ES SEED LS ES NS ES LE A SR LL ARLE TS LED LL AES LEED ES LOD AY LOD ARS LE ATO LED LES Nr SLL SES LRG ODA LG ES: I LLG ALLL LAID LILLE LD SESS OAD I, SANA Ny 


DATA BLOCK E 


ESCRIPTION: Material Properties 


Format and Example 


| ¢ | RHO | SIGMIN | SIGMAX | KEULER | FORMAT | 
| | | | | | SF10.0 | 
{ [ | 


ee 
{1000000.).. .1 | -15000.| 20000. | 4. | 
| Se eee ae 





E-Young's Modulus 

RHO~ material density 

SIGMIN-minimum allowable stress 

SIGM AX-maximum allowable stress 

KEULER~Euler buckling coefficient 
for bar element 


wm & WwW WN = 


NOTE: Number of cards read=NMT 


i 
| 
| 
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ae ae ae ee ee oo eS OE GE eg cs ES ei es a ee ee a om ee ee ee a ee ee ee 


Sa ee ees | 


| 


piensa UO ey as ey es Cy i RN CR i cil, > CS pee See pe YA ge i Se fe lem aera A a 5 I gy MU 


— 


DATA BLOCK F 


=e oS aE MD EE ee SS = 


DESCRIPTION: Bar Element Information 


Format and Example 


eae PS 
| ELNO | NODE1] NODE2 | MAT COD| NDSG if A |FORMAT | 
PL free 
| a Se SEEN a ee eee 
ah 2 ~ ee 
ae | a a 
ae acoenne> 1 7 #| 2 
| oP | | 2 | 1. 
Oe Ce Se Se | 
——| we, o> 1 + tT ia 
i to 4 | 2 1 } 1 te 
et ae ee ee, 
FIELD CONTENTS 
1 ELNO-~element number 
2 NODE1-global number associated with 
element node 
5 NODE2-~global number associated with 
element node 
4 MATCOD-material type of this element 
5 NDSG-design variable number associated 


with this element 
6 A-member cross-sectional area 


NOTE: Number of cards read=NE 
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| 
| 


ee ee a nt 


. 





Speers | iz {| 








‘stl 


ONTENTS 


JN-joint number 


IX 
sey 
IZ 


Rev ,z Constraint identifier. 


zero the corresponding degree of freedom 


is constrained) 


NOTE: Number of cards read=NCJ 
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| FORMAT | 
1 yrs 
|! Sa 


(Lf non- 


ee 


I SD ee EE ei ee ce ee Se ee ee se SE ce ee Se ee eee ee ee ee ce cee ee OS ee Dc ee ee ee SS ee a el ss ee Ge eis i) ee 





i ice cy ey a ee ES ce Oe, I LE EE LN LOO LN LS ENS LLG LLY A LA A DY ES A ED A ED LS LD A ES TS I ES A TID ce SD COD eas SE 


az 


g 





DATA BLOCK 4 


OMIT this card 


cy 


a 


| vat4o. | = 14140. 7 —— - | 





2 [tw wo.P -t4t0.f 0. TT 


Ks 
t 
tx 
t 
7] 


- WN = 





nn con 





_s 


if NLC=0 was read in block B. 


Loading Information 


115, , 3F 10. 0) 


| 
—— aa 








CONTENT 


NLJ-number of loaded joints 
for this load condition 
JN~- joint number 
FX 
PYoerorces in ~hne X,Y,2Z directions 
FZ 


Ue 
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1 


I aie cI RE Gs gg I CN ies OI Se ee A Ge, I a I ce em a ccm gp imi MI My CE ce a li a ee 


- 


IN MASS FORMAT 
| | }IS,F10.0 | 
\ 
ae aa ee 
p 1 500. | | 
can os ie ae pone hi 
FIELD CONTENTS 


OMIT this block if NFMASS=0 was read in block B. 


ESCRIPTION: Lumped Mass Information 


Format and Example 


JN-joint number 
Z MASS -concentrated mass at joint (JN) 


mass is in force units 


Le a ee ee | 
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| 


ee AU es OL cece VE sees il eel eninge OS) semi, We eee AMEE cia En emia AU) ERS ae ie AED ES A A ce ea ARN A AD ieee ce ND OR: eee eee AN GP I es, ee ee, A OUR es) CEE ees, SEE Go . ee ES wee OI ely OE ee es 


> | 
3 
to 
log 
ir 
1O 
a 
Cy 


OMIT this block if NDVAR1=0 


DESCRIPTION: Design Variable Information 
(AREA Variables) 


Format and Example 





at # - | FORMAT | 
| ( 8F10.0 | 
| = ae eee. 2 Att 
Tt. oT : ——— 
{ veer eh meline ee 
[ 01 01 | sme XAL(Z) 
ee a 
[ 10. 10. a - XAU (I) 
[ea ee 
FIELD CONTENTS 

1 XAtinitial value of area design variables 


XAL~lower bounds on area design variables 


XAU~-upper bounds on area design variables 


NOTE: read one value of XA, XAL, XAU for each 
independent area variable defined in Block D 


X 


Number of cards read =as required 


— 
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AB Sy mR «iy I RU A A Ry AAS I it SS! Ra I) AS ET OR Ay COE I I EE RN RI 5 a RS A i aM I LY OPM iy MONS I iain cs ee a a 


| 
| 


| 
| 


OMIT this block if NDVAR2=0 


DESCRIPTION: Design Variable Information 
(COORDINATE Variables) 


Format and Example 


eee nes (ee ll 
XC (1) Gel) Wes 3. « Capel ani | 


poe 








| 
; 


al 
a 


XC-initial value of coord. design variables 
XCL-~lower bounds on coord. design variables 
XCU-upper bounds on coord. design variables 


NOTE: read one value of XC, XCL, XCU for 
each independent coordinate variable 
defined in Block D 


Number of cards read =as required. 
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o-- onl 


ct 


| 





DATA BLOCK L 


OMIT this bleck if NDSPLC=0 in Block A 


DESCRIPTION: Joint Displacement Constraint Information 


Format and Example 


= | af | 





| 
| 
| 
{ 
| 
| |315, 2F10. 
| ‘ae ee 
oe ee 

rnin we |] -.02) .02 | 

—— a Sa | { 2 ' e 
{ orem ranean eaten 
| | ie | 2 1 | -.02] .02 | 
{| FIELD CONT ENTS 
| 1 JN-joint number 
2s DIR-direction 1=X 2=Y 3=2 
| O=limit on resultant displ. 
| LO-loag conda tion 
| 4 BL-lower bound on displacement 
(Tf DIR=0 read O here) 
| > BU-upper bound on displacement 
7 Number of, cards read= NDSPLC 
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—ia ee SD cee) eS ce ees SS 9 © ee, ee eee eee ce i OG ee ee ee ce ee es a ee ee 





ei Mg ETance ei mes TE a a cc ce i rl a as I ah ye i le...) es ee 
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OMIT this block if NFREQ=0 was read in Block A 


DESCRIPTION: Frequency Constraint Information 


Format and Example 


FC1 C2 | FC3 (FCS wees ees PCM | }roRmat | 


Lb fb pres } 





— Berea aay ree re SN 
20. 
eee eeemeaee 
FIELD CONT ENTS 
1 FC1- lower bound on first natural frequency 


constraint in CPS. (cycles per second) 


N FCN- lower bound on NFREQ-th natural 
frequency constraint in CPS. 


(cycles per second) 


NOTE: OMIT this block if NFREQ=0 was read in data 


block B. 


Number of cards read= as required 
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ioe ecal EXAMPLES 


A. INTRODUCTION 

Design of planar trusses and 3-dimensional space towers 
are presented here and the corresponding numerical results 
are summarized to demonstrate the use of the code. 

In the examples given here, the design variables are 
member cross-sectional areas and joint coordinates. It 
should be emphasized that in practical design, the reciprocal 
of the member areas is usually a better choice for the design 
variables. However, the purpose here is to knowingly create 
difficult optimization problems, thus the choice of variables. 


The examples begin with the 3-bar truss. 


B. CASE 1: 3-=BAR PLANAR TRUSS 

The simple 3-bar planar truss, as shown in Figure 6.1 has 
been previously used for the user guide example. This struc- 
ture was designed for optimum geometry subject to a set of 
two load conditions, buckling constraints, displacement con- 
straints, and a lower bound on the system's first natural 
frequency. The allowable stresses specified are 

-15000. < ao, < 20000. psi ae 193 

The member areas were linked in the following groups: A1=A3; 


Le 
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Woe ase a 
Only the stress constraint was imposed for this case. 
Final design information is given in Table VI. The number of 
analysis for this design was 38. 
Ze. S@ase ib 
This case includes the previous constraint plus con- 
straints on buckling (KEULER=4.0), displacement (~.02 -.02 in.), 
and first natural frequency Cw =20.0). Results are given in 
Table VII. The number of analysis required for this design 


Mets 16 with 1.73 seconds of CPU time. 





aN 


14140 lbs. 14140 lbs. 


Figure 6.1 3-BAR TRUSS 


Ly 





C. CASE 2: 18-BAR PLANAR TRUSS 

Peeemierliver truss, as shOwn in Figure 6.2, has been pre- 
viously used as a standard test case for peter ural design 
[Ref. 2]. The structure was analyzed for a Single set of load 
conditions with allowable stresses being 


-20000. < o. < 20000. psi i=1,18 


Young's modulus was taken as v0 pSi with a material density 


few. l W./cu in. 





Figure 6.2 18-BAR TRUSS 


US 


i! 





The independent coordinate variables were taken as X3, Y3, 
X5, Y5, X7, Y7, X9, Y9. The member areas were linked as 
follows: Al=A4=A8=A12=A16; A2=AG6=ALO=AL4=A18; A3Z3=A7=ALL=A15; 
A5=A9=A13=A17. There are a total of four independent area 
variables and eight coordinate variables. 

1. Case 2a 

This design analysis includes stress, buckling (KEULER= 

4.0), displacement (-10.0 to 10.0 in.), and a first fundamental 
frequency of 3 Hz. Additionally, a nonstructural mass of W=5000 
lbs. was attached to node 1. The number of analyses for this 
design was 62. Results are presented in Table VIII. 26.72 


Seeonds of CPU time were used for this calculation. 


Breeeeeenor 3: 25=BAR SPACE TOWER 





200 @// ne 
, Se —— 200 3 


Reese Gas 2o-BAR SPACE TOWER 
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The 25-bar space tower shown in Figure 6.3 was designed 
for two independent load conditions given in Table X. The 
allowable stresses were specified as follows 
-40000. < o, < 40000. psi e125 
Young's modulus was selected as 9 | psi with a material density 
p=.1 lb/cu. in. Members are assumed tubular with a nominal 
@eameter to thickness ratio of D/t=100 giving KEULER=39.274. 
Symmetry was imposed in both the x-z and y-z planes. Non- 
structural masses of W=500 lbs. were attached at nodes 1 and 
2. Coordinate variables were X4, Y4¥, Z4, X8, and Y8 with the 
remaining coordinates linked to maintain symmetry. Area vari- 
ables were linked in the following manner: Al; A2=A3=A4=A5; 
A6=A7=A8=A9; Al0=Al1l; Al2=A13; Al4=A15=A16=A17; A18=A19=A20=A21; 
A22=A23=A24=A25. 
1. Case 3a 
stress, displacement (-0.35 to 0.35 in.), Euler buck- 
ling and first natural frequency limited to a value eee 16 
HZ. were imposed for this case. Final design information is 
given in Table XI. The number of analyses for this design was 


62 with 27.99 seconds of CPU time used. 


E. CASE 4: 234-BAR SPACE TOWER 
The inital layout of the tower is shown in Figure 6.4 
stress limits were as follows 
-15000. <a, < 20000. psi 121,234 
Young's modulus was chosen as steel 3.x10/ pDS1l with a material 


Sowemeyeot aluminum o =.1 lbs/cu. in. 
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1. Case 4a 
Stress constraints as well as constraints on Euler 

buckling were imposed. The resulting problem had 56 area 
design variables and 42 coordinate variables. This problem 
meguared Ji structural analyses using 72 minutes of CPU time. 
The objective function was evaluated 89 times and gradients 
were calculated 30 times. Results are presented in Tables XII 
through XIV. Although the weight of the structure was only 
reduced by 2%, it should be noted that constraints were ini- 
tially violated and the optimizer overcame the constraint 
violations. It is believed the optimum that was reached is 
not the true optimum because of the extreme non-linearity of 


the problem. 
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Wapile. Vi 


cep enwiso. — DEsmeletNPORMATION (STRESS) 











BODY 

LOAD JOINT LOADS 

eexD. xX ne Z 
i YZ Ceetiereteo =—O.,1414E+05 0.0 
2 2 Opt ees =—O5TFteE+05 0.0 


fee MASSES 
NO. JOINT eS 


i 2 500 lbs. 


AREA (sq. in.) 


Peele, R er Pie FINAL 
Al=A3 ieee ,16602E200 
A2 iL ee ~4833E+00 


Coon Pa Aihas Cans ) 
JOINT INITIAL FINAL 


a 10 Aeros eveau al 


FINAL WEIGHT FOR AREA AND COORDINATES 


WEIGHT = 2.6349 (l1bs.) 


20 





BODY 


JOINT 
NUMBER 


Z 


2 


MEMBER 
Al=A3 


A2 


JCal 


Fable VII 


3-BAR TRUSS. DESIGN INFORMATION 
(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 


JOINT LOADING INFORMATION 


SAME AS TABLE VI 


JOINT DISPLACEMENT INFORMATION 


Dibk. LOAD LOWER Oe agean 
COND. BOUND BOUND 
il at -.02 OZ 
2 Je ~.02 OZ 


Piene SGie- Tl.) 


INITIAL FINAL 
i .30 GPSS) She ionae 21816) 
io 0.1000 E-0O1 


COORDINATES Cin.) 
INITIAL PINAL 


HO. 0 BUS ISHS evoisiasi0l ah 


PINAL WEIGHT = 2.8381 (l1bs.) 


CRITICAL CONSTRAINT - FREQUENCY 


on 





iu DY 


LOAD COND. 


NO. 


Basie. VITT 


bPooean MnUsas. GOAD CONDITIONS 


LOADM@ONDITIONS Clibs.) 


Oh 
A 
aL Oro -0 
Z = 0 -0 
~ Or -0 
6 0.0 -0 
8 Oea0 -0 


PExXEDetASso INFORMATION 


SHON a 


AL 


a2 


LOADS 

y4 Z 
~eOUO0B+0 5 Bia 
Se OES Cs oo 
{ UCN elONs Oe. 0 
e000 E05 07-0 
ZOO tO 5 O20 


MASS 


SOO Leases) 





Table IX 


18-BAR TRUSS. DESIGN INFORMATION 
(STRESS, DISPLACEMENT, EULER BUCKLING, FREQUENCY) 


AREA (sq. in.) 


MEiok R INITIAL FINAL 
A1L=A4=A8=A12=Al16 nO iO Sa520 
A2=A6=A10=A14=A18 150) .0) Sorae 7 
A3=A7=A11=A15 L050 leche 
A4=A9=A13=A17 mer 0 CUO 


COORDINATES (in.) 


JOINT INITIAL FINAL 
Xx Y X 
5 15019) 0% 0 BOS. 
é To0 0. 628 
Hl 200 Oe 384 
S 250 Oy. 148 


FINAL WEIGHT = 13,940 (l1bs.) 


Serer EeeeNSTRAINT - FREQUENCY 


5) 6 





BODY 


EBOAD COND. 


NO. 


Table X 


Z2= oe oeioo. 9 LOAD CONDITIONS 


LOAD CONDETIONS (1bs.) 


IONE LOADS 

X Y 
1 O..0 2000 
2 020 -2000 
1 SECO 0) 10000 
Z 0 10000 
3 900 0 
6 200 0 


FEXED MASS INFORMATION 
eo Cun 
zi 


2 
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MASS 
a0 tree bs.) 


SU Glace») 





Table XI : 


ao a nviooeiotGN ENFORMATION 
Senso. DisPuaeribihesEUGER BUCKLING, FREQUENCY) 


AREA (sq. sligkee) 


MEMBER INITIAL FINAL 
Po ae 1.46 
A2=A3=A4=A5 ZO 0.74 
A6=A7=A8=A9 Lae 03:8 
Al0=Al1l 22.0 aes 
A1l2=Al13 ZO reel? 
A1l4=A15=A16=Al7 Zee 0.74 
A18=A19=A20=A21 220 Ove 
A22=A2 3=A24=A25 2.0 Oe. 


COORDINATES Cans) 


CONT I BERG ess, FINAL 
X xe Z X ne 

+ cS Spf ae 100 So4 46.9 

8 LOCO L002 6 O20 1 Syele: ysis 


Pei WElGHT = 267 (lbs.) 


Serica CONSTRAINT - FREQUENCY 


aS 


L33o.3 


dh 





BODY 

LOAD COND. 
ih 
iL 
2 
2 
é 
3 
m 
m 
5 
5 

NO. 

A 

2 


Table XII 


234-BAR SPACE TOWER. LOADING INFORMATION 


ECAP RGONDERERONS Clbs .) 


Je EN Ge LOADS 
xX Be 
4 9 6000. -20000. 
20 6000. -20000. 
49 ou Cos. -20000. 
20 -o000. -20000. 
49 5) 0 OO -20000. 
50 3000. -10000. 
49 3000. -10000. 
50 G19) 0) 2 -10000. 
49 -3000. BEC) 08) Oke 
50 soil) Wie ILG(OG OFA 


Peo DeMooo lL NFORMAT TION 
JOINT 
49 


S10 


96 


MASS 
ZOCr Clinch 


2007 Cilbse 


i! 





Table XIIil 


234-BAR SPACE TOWER. DESIGN INFORMATION 


AREA (sq. in.) 
MEMBER 


A1L,A2,A3,A4 
AS,A6,A7,A8 

#99A10,A11,A12 

Ak3,A14,A15 ,A16 

A17,A18,A19,A20 

Pam A22.A23,A24 

A25 ,A26,A27,A28 

A29 ,A30,A31,A32 

A33,A34,A35 ,A36 

A37,A38,A39 ,A40 

Pie, At? AN 3, A44 

A4¥S ,A46 ,A4U7 AUB 

A49 ASO ,A51,A52 

A53,A54,A55,A56 

AS7,A58,A59,A60 

A61,A62,A63,A64 

A65 ,A66 ,A67,A68 

A69,A70,A71,A72 

A73,A74,A75,A76 

Pog A78,A79 ,A80 

A81,A83 

A82,A84 

A85 ,A87 

A86 ,A88 

A89 ,A90,A91,A92,A93,A94 ,A95 ,A96 
A97,A98,A99,A100 ,A101,A102,A103,A104 
A150 ,A106,A107,A108,A109 ,A110 ,A111,A112 
Pups AWS A115 A116 ,A117,A118,A119 ,A120 
Pan ed? .8173,A124 A125 A126 ,Al127.A128 
A129 ,A130 ,A131,A132,A133,A134,A135,A136 
Pape?) A138, A139.A140,A141 A142, A143,Al44 
A145 ,AL46 ,A147,A148,A149 ,A150,A151,A152 
A153,A154 ,A155 ,A156,A157,A158,A159,A160 
A161,A162,A163,A164,A165 ,Al166 ,A167,A168 
A169 ,A170,A171,A172,A173,A174,A175 ,Al176 
Aly 7.Av 78,4179 .A180 

Pale Aug A183,A184 

A185 ,Al86 ,A187,A188 

A189 ,A190 ,A191,A192 


er 


LO ESB 


) (65) 1 1S) (SS) 1S eS) 1S SS SSS 2 ee eer SS 1 ef e Fe Fee ees Be eke te eke Re eens 





Meiete XITE (Cacont'd) 


AREA (sq. in.) 


MEMBER ia) eee A 
A193,A194,A195 ,A196 25.0 
A197,A198,A199,A200 25.0 
A201,A202 25-0 
A203,A204 C50 
A205 ,A206,A207,A208 25.0 
M209 .A210,A211,A212 25.0 
PO13,A214 25.0 
A215 ,A216 25.0 
A217,A218 25.0 
Pes), A220 25.0 
Poo A222 25.0 
W223 A224 25.0 
A225 ,A226 a 0 
FOOT. A228 25.0 
A229 ,A230 25.0 
A231,A232 25.0 
A233,A234 25.0 


COORDINATES (in.) 


JOINT INITIAL FINAL 
X Me Z X x 
1 120 0 We AZ U. 25 0 
S a 120 ass WAC Oy D2 Oe 4 
~ FOZ 240 ey. ORs. Z 2 NAL A Z 
3 hc) 360 a3 Orel Soyo 
1ef 84 4 80 84 84.0 ie eG 
ZAM, Is 600 1S) Theo 604.5 
Zo 66 720 66 Broie 724.2 
Zo 2) 840 a7, oS aie 842.6 
ore) 48 cee 48 47.8 960.4 
Initial weight = 84,524 lbs. Final weight = 83,142 lbs. 


*Final weight (AREA VARBIABLES ONLY) = 28,900 lbs. 
* This clearly indicates optimum not reached with previous case. 


S8) 


18) 
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ieee XI TI (cont'd) 


worNT 


S)/ 


4} 
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49 


ak 


Spe) 


COORDINATES (in.) 


ENITIAL 
x 


1080 


*“NOT A DESIGN VARIABLE 
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60 
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Low 


Bs) 


FINAL 
a4 


O77 .4 


os ee 


hele 





i eeOnNc DUS AND KECOMMENDATIONS 


eee CONCLUSIONS 

The finite element code was presented coupled to an 
Spummazer for truss analysis and design. Trusses were 
designed for minimum weight with multiple load conditions 
considered. 

The displacement method for static analysis and the sub- 
Space iteration method for eigenvalues were applied. 

several examples were considered. In every case the code 
worked as an analysis tool, and significant weight reductions 
were obtained with the coupled optimizer CONMIN. Run times 
and results compared to the same test cases with other codes 


indiate that the code is competitive as a design tool. 


B. RECOMMENDATIONS 
The following recommendations may be of value for follow 
on work: 


1. The code should be modified so the user can access 
stresses directly. 


2. The code should be extended to include other elements 
such as frames and plates. 


3. An out of core equation solver should be added. 
4. The method of gradient calculation should be dependent 


on specific gradients required reference 3 and reference 
4, 
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Gradients of frequency constraints would benefit from 
a more efficient algorithm reference 5. 


The need for a large scale public structural optimiza- 
Enron Code still exists. 
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APPENDIX A 


DATA FILES 


A. INTRODUCTION 
This appendix contains the data files used to create the 
test cases in Chapter VI. Additionally the data file for the 


user's guide in complete form is presented. 





fap le xi V 


DATA FILE 3-BAR TRUSS 


CStRes 5) 


3="GrRhe TRUSS Chest CASE A=-STRESS ONLY) 
& 4 1 2 Q 0 


3 0 
0 0 1 
0 386% 0.0 
l -10. Oe 0. i 
2 06. -10. Oe 0 
5 O6« O06 0. 0 
4 10. O-« 0. Z 
10000000. el ~15000. 20000. 4e 
sh i 2 1 1 le 
2 S Z lL 2 le 
3 % 2 Ll 1 le 
1 l i 1 - 
2 Q 0 1 
3 l l l 
: ik 1 l 
‘ 14140. -14140. O6« 
2 -14140. -14140. OQ. 
, 2 500-6 
Paepe eO1 
10. 10. 
10. 
1.0 
206 
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Table XV 


DATA FILE 3-BAR TRUSS 


(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 


— ee 


3 4 
2 2 L 
0 38664 0.0 
L -10. Oe QO. 1 ) 0 -l. 
2 QO. -10. Or 0 0 6) QO 
3 O~« O. O-« 0 0 9) OQ. 
+ 10. O~« 0. 1 8) 9) le 
l 1 2 1 l le 
Za 2 2 L 2 le 
3 & 2 L 1 le 
l 1l L L 
2 0 @) 1 
3 l lL 1 
: 1 1 1 
: 14140. -14140. 0. 
2 -14140. 7-14140. 0. 
2 500. 

le le 

QL Ql 

10. 10. 

10. 

1.0 

206 
2 1 1 002 2 
2 2 lL -.02 02 

20- 
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XVI 


Table 


Vee Dees, Rk RUSS 


(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 
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Table XVII 


DATA FILE 25-BAR SPACE TOWER 


(STRESS, DISPLACEMENT, BUCKLING, FREQUENCY) 


NOW 
a Oo 
Oo 
ce) 
e 
° eeete Oo 
QO aAHtAtHO 
=@ 
°° ee oe 
© Confe~d © 0 ~de-4 700 
QO tmtt Ital ft NOW 
(@] 
e 
(eo) 
oe om 
e ee ° © e 
oO ®@ © Or-ted © Or-4 “OO 
uw CQ tedet | f mtet NOW 
fod t ce) 
LY & r 
~ oO 
es) ‘ 
Zz OoOMAMNMdOOO9O 
i 
—J wt onal e ° 
» 4 toed Oo7Oo 98 
a ) Ca | 20 O0DHO 
ae | OQOKRINAINAIN WNW ¢« Nowe en 
co oO Oo 
~ om e 
° Oo 
a. 
MNO DO At Het SEES 
m 
Q ; 0:0 0 
~~ ® 
Ww? 4 2000H0 
vt vo) NOINA oth 
us oO Oo 
oc eee © @ t+ °° ee @ 
i= Co000090o0 oo oo Oo WIL UY 
v”) oOn0000 © eee eee @eeeveeeveeet®teevsensvse? ee 8 06 8 @ Oo VWOee = MAM 
ond NAMA HHODOO ANAANNANANANAANNNNNAN ANNAN 9oo0o 909090 eee 
y i ne 
Vv) e 
<< ° O20 Oo 
oS ] © eO00HO 
NIA OD ANANGISM MOM & FINN OO OOP hm 20000 © NOue eh 
b= in fo] >] 
WY) WV et © ¢ © 8G °° eg TALALALY 
uw OQeem eO000 ¢ «*e OO Oo MAM Mg 
t ee eM MO-O000 1 Q0 OO ee eeeoe 
ww A HOOOMM I rel paw poe pew ee ee re eee ee ee ee ee ee ee ee ae ee ee ee eee) od etitiend a2 teed ee 
Vv) ; O0 ct wa WN ) 
Vv} NIN Oer & 
baie = | 20OrH+O 
con ONES COUN DO FIN OD 0 1909.0 OF FIN ON. O09 FIN i NOWM eiNete-ietea 
est + wn ° e ~ ° oO 
e NM AMNINNO e © oO ee ee ® 
tw ON e& © © eO000 Oo 9oO0090 Oo 
es) DO ee MRE M HOO ee O 90000 
i ia mala Pim a hed naam an DNAS ead mae Soe aaetet OO AOMNININ mt ANIA 
a) ] oO pea] 6 ° 
4 © e © | 
oc eOONHO 
LINO NOM OTD MVE NON ONO AINE HH ON DP IAN FAAS NORD PONANE ANION GING tery 
ON =19 ond end rnd 4 0h eh A Ht ANI CIO ONIOIO °o 
e 
1 = 3 
N 


16.0 


wee 
. — 


66 





Table XVIII 


DATA FILE 234-BAR SPACE TOWER 
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234-BAR SPACE TOWER DESIGN (TEST CASE STRESS »BUCK. »yDISPL.» FREQ.) 
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PerewoLA 3 


PROGRAM ORGANIZATION 


Pee beockiPTiONn 

The program organization is laid out in the following 
flow charts. The main driver program is arranged to call a 
pueariver, XAMSADT, and the optimizer of the user's choice. 
All changes for replacing an optimizer occur in MSADT. This 
allows for easy testing of several optimizers on the same 
problem. 

XMSADT may be called from the main for input, analysis, 
and output. Printed output may vary as the user requires. 
A complete listing of all subroutines and their functions is 


given in Table 
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